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Abstract 

A general form factor formula for the scaling Z(N)-Ismg model is constructed. 
Exact expressions of all matrix elements are obtained for several local operators. 
In addition, the commutation rules for order, disorder parameters and para-Fermi 
fields are derived. Because of the unusual statistics of the fields, the quantum field 
theory seems to be not related to any classical Lagrangian or field equation. 
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1 Introduction 



The 'form factor program' is part of the so-called 'bootstrap program' for 
integrable quantum field theories in 1+1-dimensions. This program classifies 
integrable quantum field theoretic models and in addition it provides their 
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explicit exact solutions in terms of all Wightman functions. This means, in 
particular, that we do not quantize a classical field theory. In fact the quantum 
field theory considered in this paper is not related (at least to our knowledge) 
to any classical Lagrangian or field equations of massive particles. The reason 
for this seems to be the unusual anyonic statistics of the fields, turning this 
form factor investigations even more fascinating. The bootstrap program con- 
sists of three main steps: First the S-matrix is calculated by means of general 
properties as unitarity and crossing, the Yang-Baxter equations and the addi- 
tional assumption of 'maximal analyticity'. Second, matrix elements of local 
operators 

out (p' m ,...,p' 1 \0(x)\p 1 ,...,p n ) m 

are calculated using the 2-particle S-matrix as an input. As a third step the 
Wightman functions can be obtained by inserting a complete set of interme- 
diate states. 

The generalized form factors [1] are defined by the vacuum - n-particle matrix 
elements 

( | 0(x) | Pl ,..., Pn )Z... an = e-^+-+^ C.^, . . . ,0 B ) 

where the ctj denote the type (charge) and the 9{ are the rapidities of the par- 
ticles (pi = Mj(cosh#j, sinh#j)). This definition is meant for 9± > ■ ■ ■ > n , in 
the other sectors of the variables the function F®(ff) = F® (0i, . . . , n ) is 
given by analytic continuation with respect to the 0^. General matrix elements 
are obtained from F®(0) by crossing which means in particular the analytic 
continuation 0, — > Oi ± in. Using general LSZ assumptions and maximal ana- 
lyticity in [2] the following properties of form factors have been derived: 4 

(0) The form factor function F®(9) is meromorphic with respect to all vari- 
ables 0i, . . . , n . 

(1) It satisfies Watson's equations 

F °.a iaj ...(- • • > 0j> • • • ) = F °. ajai ...(- ■ ■ > Qji hi ■ ■ ■ ) S a i a j (0ij)- 

(ii) The crossing relation means for the connected part (see e.g. [4]) of the 
matrix element 

| o(o) | P2 , . . . , Pn y^T n = F° ia2 ...* n (0i + ™, o 2 ,..., n ) 

= F^ anai (9 2 ,...,9 n ,9 1 -z7r) 

where a°(a) is the statistics factor of the operator O with respect to the 
particle a. 

4 The formulae have been proposed in [3] as a generalization of formulae in [1]. The 
formulae are written here for the case of no backward scattering, for the general 
case see [4]. 
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(iii) The function F®{9) has poles determined by one-particle states in each 
sub-channel. In particular, if a\ is the antiparticle of a 2 , it has the so- 
called annihilation pole at 9i 2 = in which implies the recursion formula 

Res F^e,, . . . A) = 2i c aia2 F°^e 3 , . . . A) 

d 12 =lTT — J 

x (l - o-°(a 2 )S aian (9 2n ) ■ ■ ■ S a2a[i (6 23 )) ■ 

(iv) Bound state form factors yield another recursion formula 

ResFS...^, e 2 , e 3 ,...) = V2i£.(0 ( i2) A, ■ ■ ■ ) rjj 

if z?7 is the position of the bound state pole. The bound state intertwiner 
T^g (see e.g. [5,4]) is defined by 

iResS Q/3 (0) = rfr^ 

6=IT) ' ' 

(v) Since we are dealing with relativistic quantum field theories Lorentz co- 
variance in the form 

F^(9 1 ,...,9 n ) =e a "F| , (0 1 +//,..., 8 n + ») 

holds if the local operator transforms as O — > e Sfl O where s is the "spin" 
of O. 

Note that consistency of (ii), (iii) and (v) imply a relation of spin and statistics 
a°(a) = e~ 2ms and also a°(a) = 1/a (a) where a is the anti-particle of a, 
which has the same charge as O. All solutions of the form factor equations (i) 
- (v) should provide the matrix elements of all fields in an integrable quantum 
filed theory with a given S-matrix. 

Generalized form factors are of the form [1] 

II HO*), (* y = *i-0;) (1) 

l<i<j<n 

where F(6) is the 'minimal' form factor function. It is the solution of Watsons 
equation [6] and the crossing relation for n = 2 

F(8) = F(-0)S(0) 
F(m -6) = F(in + 6) 

with no poles and zeros in the physical strip < Im 9 < n (and a simple zero 
at 9 — 0). In [4] a general integral representation for the K-function K®(9_) in 
terms of the off-shell Bethe Ansatz [7,8] has been presented, which transforms 
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the complicated equations (i)-(v) for the form factors to simple ones for the 
p-functions (see [4] and (41) below) 

K°(6) = J2c nm [ d Zl ---[ dz m h(9,z)p°(9,z)^^9,z) . (3) 

The symbols Cg denote specific contours in the complex Zj-planes. The function 
h(9, z) is scalar and encodes only data from the scattering matrix. The function 
p°{9_, z) on the other hand depends on the explicit nature of the local operator 
O. We discuss these objects in more detail below. For the case of a diagonal 
S-matrix, as in this paper, the off-shell Bethe Ansatz vector ^a{9_, z) is trivial. 
The K- function K®(9) is an meromorphic function and has the 'physical poles' 
in < 1m 9ij < n corresponding to the form factor properties (iii) and (iv). 
It turns out that for the examples we consider in this paper there is only one 
term in the sum of (3). 

In this paper we will focus on the determination of the form factors of the 
scaling Z(iV)-Ising quantum field theory in 1+1 dimensions. An Euclidean 
field theory is obtained as the scaling limit of a classical statistical lattice 
model in 2-dimensions given by the partition function 



^exp 



2wi /N 



M \ (ij) 

as a generalization of the Ising model. It was conjectured by Koberle and 
Swieca [9] that there exists a Z(7V)-invariant interaction E(a i: aj) such that 
the resulting massive quantum field theory is integrable. In particular for 
N = 2 the scaling Z(2)-Ising model is the well investigated model [10,11,12,13] 
which is equivalent to a massive free Dirac field theory. In this paper we in- 
vestigate the general Z(iV)-model. It has also been discussed as a deformation 
[14,15] of a conformal Z(N) para-fermi field theory [16]. The Z(N)-lsmg model 
in the scaling limit possesses N — 1 types of particles: a = 1, . . . , N — 1 of 
charge a, mass M a = Msin7r-^ and a = N — a is the antiparticle of a. The 
n-particle S-matrix factorizes in terms of two-particle ones since the model 
is integrable. The two-particle S-matrix for the Z(iV)-Ising model has been 
proposed by Koberle and Swieca [9] . The scattering of two particles of type 1 
is 

_ sinhl(0 + 2*) 
SW -si„h (4) 
This S-matrix is consistent with the picture that the bound state of iV — 1 
particles of type 1 is the anti-particle of I. This will be essential also for 
the construction of form factors below. We construct generalized form factors 
of an operator 0(x) and n particles of type 1 and for simplicity we write 
— F\ i($)- Note that all further matrix elements with different particle 
states of the field operator 0{x) are obtained by the crossing formula (ii) and 
the bound state formula (iv). As an application of this form factor approach 
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we compute the commutation relations of fields. In particular, we consider 
the fields iPqq{x) , (Q, Q — 0, . . . , N — 1) with charge Q and 'dual charge' Q. 
There are in particular the order parameters ctq{x) = i/jqq(x), the disorder 
parameters (J>q(x) = ^oq( x ) and the para-fermi fields i/)q(x) = i^qq{x). We 
show that they satisfy the space like commutation rules: 

<j q( x ) <j Q' (y) = VQ'(y)vQ(x) 

^Q( x )^Q>(y) = ^Q'(y)^Q( x ) 

Hq(x)*q(v) = a Q {y)^{x)e e ^-y^^/ N (5) 

M*)Mv) = My)M x )e e(xl ~ yl)2mQ2/N • 

It turns out that the nontrivial statistics factors a°(a) in the form factor 
equations (ii) and (iii) lead the nontrivial order- disorder and the anyonic 
statistics of the fields. 

The form factor bootstrap program has been applied in [13] to the Z(2)- 
model. Form factors for the Z(3)-model were investigated by one of the present 
authors in [17]. There the form factors of the order parameter o\ were proposed 
up to four particles. Kirilov and Smirnov [18] proposed all form factors of the 
Z(3)-model in terms of determinants. Some related work can be found in 
[19]. For general N form factors for charge-less states (n particles of type 1 
and n particles of type N — 1) were calculated in [20]. In the present paper 
we present for the scaling Z(7V)-Ising model integral representations for all 
matrix elements of several field operators. 

Recently, there has been a renewed interest in the form factors program in 
connection to condensed matter physics [21,22,23] and atomic physics [24]. In 
particular, applications to Mott insulators and carbon nanotubes as well as 
in the field of Bose-Einstein condensates of ultra-cold atomic gases have been 
discussed and in some instances correlation functions have been computed. 

The paper is organized as follows: In Section 2 we construct the general form 
factor formula for the simplest N = 2 case, which corresponds to the well- 
known scaling Z{2) - Ising model, and show that the known results can be 
reproduced by our new general approach. In Section 3 we construct the gen- 
eral form factor formula for the Z(3) case, which is more complex due to 
the presence of bound states, and discuss several explicit examples. We ex- 
tend these results in Section 4, where the general form factors for Z(N) are 
constructed and discussed in detail. Section 5 contains the derivation of the 
commutation rules of the fields and some applications of this formalism are 
presented. Some results of the present article have been published previously 
[25] without proves. 
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2 Z(2)-form factors 



To make our method more transparent and with the hope that our construc- 
tion will help to calculate form factors for all primary and descendant fields, we 
start with the simplest case N = 2, which corresponds to the well known Ising 
model in the scaling limit. This model, already investigated in [10,11,13,20], 
is equivalent to a massive free Dirac field theory. The model possesses one 
particle with mass M and the 2-particle S-matrix is S = —1. In [13] the form 
factor approach has been applied to this case with the result for the order 
parameter field o~(x) 

FZ(9) = (2i)^ 2 J] tanh l% ( 6 ) 

l<i<j<n 

for n odd. It is easy to see that this expression satisfies the form factor equa- 
tions (i) - (iii) with statistics factor a a = 1. For the Z(2) case in this section 
we skip the proof that the functions obtained by our general formula satisfy 
the form factor equations (i) - (v), the reader may easily reduce the proofs for 
the Z(3) and the general Z(N) case of the following sections to this simple 
one. Rather, we present the results for several fields, in particular, we will 
show that our general formula reproduces the known results. 



2. 1 The general formula for n-particle form factors 



We propose the n-particle form factors of an operator 0(x) as given by (1) 
with the minimal form factor function 

F(0) = sinhi0 (7) 

which is the minimal solution of Watson's equations and crossing (2) for S = 
— 1. The K-function is given by our general formula (3) where the Bethe Ansatz 
vector is trivial (because there is no backward scattering) and the sum consists 
only of one term 

K°{9) = N n I nm {d )P °). (8) 
The fundamental building blocks of form factors are 




fcM=nn^-^) n m 

i=l j=l l<i<j<m 

The h-function does not depend on the operator but only on the S-matrix of 
the model, whereas the p-function depends on the operator. Both are analytic 
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functions of 0j (i — 1, . . . , n) and Zj (j — 1, . . . , to) and are symmetric under 
9i <-> and ^ <-> z^. For all integration variables the integration contours 
Cg_ = Y^Cqi enclose clock wise oriented the points Zj = 9i(i = 1, . . . ,n). This 
means we assume that <fi(z) has a pole at z = such that i? = J c dz<j>(z — 0). 
The functions (z) and r(z) are given in terms of the minimal form factor 
function as 

1 -2i 
~ F(9)F(6 + in) ~ imhl (U) 

T ^ = -f^)=-^ Z - (12) 

The following properties of the p-functions guarantee that the form factors 
satisfy (i) - (iii) 

(i' 2 ) VnmiQ-i 1S symmetric under 6>j <-> 9j 

(ii' 2 ) o°v°JP\ + 27ri, 9 2 ,...,z)= p°J9 u 9 2 ,...,z) 

(iii' 2 ) if = ivr : p£ m (0, i)| 21=(?1 = <T°^ m (£,i)| zl=fla 

= (7 °Pn-2m-l(^3, • • • , #n, 32, • • • , Z m ) + p 

where a° is the statistics factor of the operator O with respect to the particle. 
The function p must not contribute after integration, which means in partic- 
ular that is does not depend on the Zi (in most cases p — 0). For convenience 
we have introduced the indices nm to denote the number of variables. For the 
recursion relation (iii) in addition the normalization coefficients have to satisfy 

N n = iN n _ 2 . (13) 

One may convince oneself that the form factor satisfies (i) and (ii). Not so 
trivial is the residue relation (iii), however, it follows from a simplified version 
of the proofs for the Z(3) and the general Z(N) case below. 



2.2 Examples of fields and their p-functions: 



We present the following correspondence of operators and p-functions which 
are solutions of (i'2) - (iii'2) - For the order parameter o~(x), the disorder param- 
eter (J,(x), fermi field ip(x), and the higher conserved currents Jl(x), (L e Z) 
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we propose the following p-functions and statistics parameters 



a(x) <-> V°anS3-i = 1 for n = 2m + 1 with a a = 1 

fi(x) ^ p£ m (9,z) =i m e& z ^^ e ^ forn = 2m with ^ = -1 

(14) 

^(x) <-> p#*(0,z) = e^^EO for n = 2m + 1 with ^ = -1 
(x) <-> z) = E e ±e < E e Lzi for n = 2m with o» = 1 

Note that p 7^ in (iii'2) only for jf;. The motivation of these choices will be 
more obvious when we investigate the commutation rules of fields in section 
5 and from the properties of the form factors which we now discuss in more 
detail. 



Explicit expressions of the form factors: 

Now we have to check that the proposed p-functions really provide the well 
known form factors for the order, disorder and fermi fields. In order to get 
simple expressions for these form factors, we have to calculate the integral (9) 
with (10) for each p-function separately. 

For the order parameter: Only for odd numbers of particles the form factors 
are non-zero. We calculate 

W£,l) = 2 m 1] tin forn = 2m+l. (15) 

l<i<j<n COSjl 2 Ui 3 

The proof of this formula can be found in Appendix B. The general formulae 
(1), (3), (7), and (13) with N n = i^" 1 )/ 2 then imply for n odd 

K{0) = 2 m n ^7T = ^ 1)/2 II tanh l% (16) 

l<i<j<n tUbn 2 U% 3 l<i<j<n 

which agrees with the known result (6). The proof that the integral I nm (9, 1) 
vanishes for even n and m > is simple: If m(m — 1) < n we may decompose 
for real 9i the contours Cg = Co — Co-in where Re Co goes from —00 to 00 
and Im (0j + in) < ImC < Im The shift Z{ — > Zi — in implies h(0,z) — > 
( — l) n h(9, z) such that the integrals along C and Co-ivr cancel for even n. 

For the disorder parameter: Only for even numbers of particles the form 
factors are non-zero. We calculate with p M as given in (14) 

Inm(0, P n=2 m J] tin fom = 2m . 

l<i<j<n COSjl 2 Ui 3 

The proof of this formula is analog to that in Appendix B, therefore with 
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N n = i n l 2 the form factors are for n even 

F»(6) = m n/2 II tanhi%. (17) 

l<i<j<n 

Similar as above for the order parameter one can show that the integral 
InmidiP^) vanishes for odd n and m > 0. It is also interesting to investi- 
gate the asymptotic behavior of the form factors when one of the rapidities 
goes to infinity. From the integral representation it is easy to check that 

Of course, this result follows easily from the explicit expressions (16) and 
(17). This asymptotic behavior is another motivation for our choice (14) of 
the p-function for a(x) and fi(x). 

For the fermi field: Only for n = 1 the form factors are non-zero. We 
calculate with p^ as given in (14) 

Inmilp^) = S nl e^¥ forn = 2m + 1 . 

The proof that I n m{9,p^ ) = for n = 2m + 1 odd and m > is the same as 
for the disorder parameter. Therefore with the normalization Ni = \/~M we 
obtain 

Ff {9) = (0 | m I 9) = u{9) = v^(Kf) • (18) 

The property that all form factors of the fermi field vanish except the vac- 
uum one-particle matrix element reflects the fact that ip(x) is a free field, in 
particular for the Wightman functions one easily obtains 

(0 | ^( Xl )---ij(x n ) | 0) = (0 | i/>(x 1 )---i/;(x n ) | 0) free . 



For the infinite set of conserved higher currents: Only for n = 2 the 
form factors are non-zero. We calculate with p J L as given in (14) 

> = ** ( e±e ' + e±fc ) 21 (SS^ + SSTfe ) for " = 2m 

The proof that Inm(9,p J ^) = for n = 2m > 2 is again similar as above. 
With the normalization C21 = ±iM the form factors are 

( I J±(0) I 0i, ... , 9 n ) m = T 5„ 2 2M (e ±ei + e ±f?2 ) (e iei - e Le >) — ^- 
such that as in [26] the charges Ql — J dxJ^(x) satisfy the eigenvalue equation 
(^-I> Lf?l ) \9 l} ...,9 n ) m = for L = ±l,±3,.... 
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Obviously we get the energy momentum tensor from J ±1 (x). 

The property F$ = for odd n > 1 and F£ = for even n > 2 is related to 
the fact that in the recursion relation (iii) the factor (1 — a° T7 S) is zero in 
both cases. 



3 Z(3)-form factors 



Let us now consider N = 3, which corresponds to the scaling Z(3)-Ising model. 
In this case we have two particles, 1 and 2, and the bound state of two particles 
of type 1 is the particle 2, which in turn is the antiparticle of particle 1. 
Conversely, the bound state of two particles of type 2 is the particle of type 1, 
which in turn is again the antiparticle of 2. So, our construction should take 
into account this structure of the bound states. We construct the form factors 
for particles of type 1, the others can then be obtained by the form factor 
bound state formula (iv). 



3. 1 The general formula for n-particle form factors 



In order to obtain a recursion relation where only form factors for particles of 
type 1 are involved, we have to apply the bound state relation (iv) to get the 
anti-particle and then the creation annihilation equation (iii) to obtain 

Res Res F^ U1 _ A (9 U •••) = Res. F°^{0 {12) , 9 3: . . . )V2T 

023=111 012="7 0(12)3 =Mr 

= 2^(04, ...)(l-aff[ 5(0*)) V2T (19) 

where 9(i 2 ) = \{9\ + 9 2 ) is the bound state rapidity, 77 = -|7r is the bound 

state fusion angle and T = i |Res6» = ^ (6?) | is the bound state intertwiner 
(see [5,4]). In the following we use again the short notation F® ± (9i, . . . , 9 n ) = 
F®(G) and also write the statistics factor as o-°(l) = erf. We write the form 
factors again in the form (1) where minimal form factor function 

r°° dt 2coshksinh|t / , / 9 \\ 

is the solution of Watson's equations (2) with the S-matrix (4) for N = 3. The 
constant C3 is given by (A.l) in appendix A. Similar as above we make the 
Ansatz for the K-functions 

K°(9) = N n I nmk (9,p°) (21) 
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with the fundamental building blocks of form factors 

lnmk(0,p) = — lyr / ^••• / / ^••• / ^Hi^,u)p(9,z,u) 

m\k\ Jcg_ R JCe_ R JCe_ R JCg_ R 

(22) 



n I rn k 



h(e } z,u) = n n <K*i - °i) n ^ - *o ^) 

i=l \j=l j=l ) 

x n r (%) n T ( u ij) n n 

l<i<j<m 1 < i < j < fc l<i<m l<j<k 

Again the integration contours Cg = ^2Ce i enclose the points &i such that 
R = J Ce dz(f)(z — 9). Equations (iii) and (iv), in particular (19) lead to the 
relations 

1 2 

[] <j>(9 + kirj) J] F(8 + kit]) = 1 (24) 

fc=0 k=0 

t{z)<P{z)<P{-z) = 1 , x(z)<f>(z) = 1 (25) 
as an extension of (11) and (12) for the Z(2) case. The solution for is 

<t>i z ) = — n — r_ TT7 rr (26) 

sinn -z sinn - [z + ir\) 

if the constant C3 is fixed as in (A.l). The phi-function satisfies the 'Jost 
function' property 

y = *<•>■ < 27 > 

We will now show again that by the Ansatz (21)-(23) we have transformed the 
form factor equations (i)-(v) in particular equation (19) to simple equations 
for the p-functions. 

Assumptions for the p-functions: The function p^ mfc (#, z, u) is analytic in 
all variables and satisfies: 

(i' 3 ) Pnmk{Q--> yd is symmetric under Q { <-> 9j 

(ii'a) ^Pnmk^i + 2ni,9 2 , ...,z,u)= Pn mk (9 1 ,9 2 , ...,z,u) 



(iii' 3 ) if 9 12 = 6» 23 = irj 

Pnmk(^^U)\zi=8i = V?P% m k(6,Z,u)\zi=e 2 

Ul=02 «1=03 

= ^?Pn-3m-lfc-l(^l'>«')+P 
( V ' 3 ) Pnmkii + H,Z + fi,U + H) = e Sfl p% mk {9, Z, U) 



(28) 
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where & = (# 4 , . . . , 9 n ), z! = (z 2 , . . . , z m ) and u' = (u 2 , • • • , Uk). In (ii' 3 ) and 
(iii'3) <rf is the statistics factor of the operator O with respect to the particle 
of type 1 and in (v' 3 ) s is the spin of the operator O. Again p must not 
contribute after integration (in most cases p — 0). Again one may convince 
oneself that the form factor satisfies (i) and (ii) if h(0,z) is symmetric under 
9i <-> 9j and periodic with respect to 6i — > 0i + 2m. Not so trivial is the residue 
relation (iii) which is proved in the following lemma. 

Lemma 1 The form factors F®(6) defined by (1) and (20)-(23) satisfies 
(i)-(v), in particular (19), if the p-functions satisfy (i' 3 )-(v' 5 ) of (28), the 
functions 0, r and h the relations (24), (25) and the normalization constants 
in (21) the recursion relation 

N n ( Res (f)(-0)) 2 4>{ir]) F 2 {ir,)F{2ir,) = N n _ 3 2iV2F . (29) 

\9=iri J 



Proof. The form factor equations (i), (ii), and (v) follow obviously from the 
equations for the p-functions (i' 3 ), (ii'3), and (v'3), respectively. As already 
stated, we will prove properties (iii) and (iv) together in the form of (19). 
Taking the residues Rese 23 =^ Res^^ there will be mk equal terms originating 
from pinchings for Zi and u^. We pick them from Z\ and u\ and rewrite the 
products that appear in the expression for I nmk in a convenient form, such that 
the location of the poles turns out to be separated from the general expression. 
Then, the essential calculation to be performed is 

mk f dz 2 f dz m f du 2 f du k 
Res Res I nm k{9,P = — ttt / -= / -= / -5- 

e 23 =iv 9i2-iv m\k\ Jc e , K Jc e , K Jc , K Jc e , K 

(m k 
n ^ - °i) n <K»i - 
i=2 j=2 

m k 

x n r (%) n ^iihi^-";) 

2<i<j<m 2<i<j<k i=2j=2 

3 / m k 

xii n^-»i)ii^- fl i 



n 

X 



with 



=1 \i=2 j=2 



Res Res f ^ f ^ f[ - 0,)^ ~ 0i)) 

623=iv 6i2=iri J Ce_ ti JCe_ K ~J[ 



rn k km 

x n r M n r ( M ij) x (^i - ui) n ^ - n ^ - ^PnH.,^. 

j=2 j=2 3=2 3=2 

Replacing Cg_ by Cqi where & = (64, . . . , n ) we have used r(0) = r(±iri) = 
x(0) = x{—irj) = and the fact that the z±, Mi-integrations give no n- vanishing 
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results only for (zi,ui) at (#1, 62) and (#2, O3). This is because for #12, O23 — ■> i?7 
pinching appears at = 2 , M i = #3 an d ^ = #1,^1 = #2- Defining the 
function 

m k km 

j=2 3=2 j=2 j=2 

and using the property that /(z, z) = f(z, z — irj) — 0, we calculate 

r = SsJ c ^L^fM ft (*(* - «.)««> - «.)) 

= (Res^(-e)) 2 2 (i^) f[ (^(^ijm.)) (/(92,9a) - /(9i,9 2 ) f[S(0 x )) 

We have used the symmetries of the phi-function (f)(9) = 0(6? + 27ri) = 0(— — 
irj) the Jost property (27) and 9 12 , 623 — irj — \%ix which implies that 

Res 0(#2i) Res 0(6> 32 ) = Res 0(0 13 ) Res <f>(8 21 ) = (Res0(-0) 

9 12 =ir) 923=ir) 9i2=irj 02:i=iv \8=irj 

W(0si) = 0(M0(M =0 2 (^) 



With the help of the defining equations (25) for r and x which imply 

(ft - ^)) = r (^i - - #2) = r(6 2 - z)k(Q x - z) 

= r(6 2 - z)k(z - 6 3 ) = t(8 3 - z)x(9 l - z) 
we obtain the relations for /(#2, #3) and f(9i, 62) 

3 / m k \ 

n n <K*i - *o n - *o m.^o = ^m^a^a,!/) 

i=l \j=2 3=2 j 

3 / m fc \ 

n n - *o n ^ - /(^ ^) = & ^, ^ # 3 , «') . 

i=l \j=2 j=2 ) 

Finally we obtain using the defining relation (24) for the phi-function 
Res Res I nmk (6,p°) = (Res0(-#)V 2 (^)x(^) {[ f[ (F(^))- 1 

B 2 3=ir) 9i2=iri \e=iij J i=lj=4 

X In-3m-lk-l(8!,P°) (l-tfU Sie^ 

if the p-function satisfies (iii'3). Therefore the form factor given by (1) and 
(20)-(23) satisfies (14) and (iii) if the normalization constants satisfy (29). ■ 
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3.2 Examples of fields and p-functions 



We present solutions of the equations for the p-functions (i' 3 )-(v' 3 ) of (28) 
and some explicit examples of the resulting form factors. We identify the 
fields by the properties of their matrix elements. In section 5 we show that 
the field satisfy the desired commutation rules. This motivates to propose a 
correspondence of fields <p{x) and p-functions p^ m fc(#, z, m)- 

The order parameter <jq(x): We look for a solution of (i' 3 )-(v' 3 ) with 

charge Q = 1 , 2 
spin s = 

statistics o\ Q = 1 . 

Since the fields carry the charge Q the only non-vanishing form factors with 
n particles of type 1 are the ones with n = Qmod3. We propose the corre- 
spondence of the field and the p-function: 



a Q <"»• Pnmfe = 1 with n = 3l + Q, < 



m = l + l k = l for Q = 1 
m = Z + l£; = / + lforQ = 2. 



The normalization constants N n follow from (29). 

Examples for Q = 1: The form factors of the order parameter a±(x) for one 
and four particles of type 1 are 

F^ = (0\a 1 (0)\p) 1 = N 1 I no = l 
F 1 \\M) = (0\cr 1 (0)\p 1 ,p 2 ,p 3 ,p A )T ul = N 4 h 21 (9,l) J] 

l<j<j<4 

where we calculate from our integral representation (22) 

(0, 1) = const (Y e~ 6i V e e > - l) TT - . - 1 . 

42H-, ) ^ ) ii.sinh - ir] ) S inh|(% + i77) 

This result has already been obtained in [17] where also the form factor equa- 
tion (iv) has been discussed, in particular (up to normalizations) 



Res (0 | oi(0) |pi,p 2 ,P3,P4)mi = const(0\ <n(0) |pi,p 2 ,P3 + P4 )m 

&34=27rj/3 



with 



< o I ^(0) I Pl ,P2, ps )Tn = -rrr Q .-^Aw* ,^ II 



sinh |(0i2 - 177) sinh \ (8 12 + M7) ~\ cosh |0 i3 
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where F™ n is the minimal form factor function for the S-matrix Si 2 . Further 
it has been found in [17] that 

Res (0\a 1 (0)\p 1 ,p 2 ,p 3 )™ 2 = const( | a^O) |pi +p 2 ,p 3 >22 

8 1 2=2iri/3 

with 

and the form factor equation (iii) has been checked 

Res (0|<7i(0)|pi,p2,p 3 )ii2 = const (S(6 12 )-l) ■ 



Example for Q = 2: The form factor of the order parameter o~ 2 (x) for two 
particles of type 1 is 

Fg(9) = ( | <r 2 (0) | p 1 ,p 2 )[ n 1 = N 2 I 211 (9, 1)F(0 12 ) 



where we calculate 

hii(d, 1) = const 



sinh |(0i2 - ^) sm h 1(^12 + irf) 



which agrees with the result (30) of [17]. This is to be expected because of 
charge conjugation. 

The disorder parameter /j,q(x): We look for a solution of (i' 3 )-(v'3) with 



charge Q = 



spin 







statistics <7i = uj 



where uj = e tr, ,n = 2/3. We call the number Q — 1,2 the 'dual charge' of the 
field fiq(x). Since the fields carry the charge Q = the only non- vanishing 
form factors with n particles of type 1 are the ones with n = 0mod3. We 
propose the correspondence of the field and the p-functiomQ 



(m n 
£ * - § £ ft. 
1=1 ° 1=1 

(m k n 

£ z» + £ «i - 3 £ 0, 
i=l i=l i=l 



with < 



n = 3m 
k = m 



where p 
(29). ^ 



^<2(<2 N+2)m ^ ^g a j n norm alization constants N n follow from 
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Examples for Q = 1,2: The form factors of the disorder parameter /j,q(x) 
for and 3 particles of type 1 are 

**** = < 1 ^5(0) | 0) = 1 

F"£(9) = (0\^(0)\p 1 ,p 2 ,p 3 )\ n 11 = N 3 I 3U (9,p^) J] FuiOij). 

i<«<i<3 

We calculate from our integral representation (22) 

n 

hu(0,p^) = const e <=i ]T e ±9 < ]J . t , , x 

where the upper sign is for Q — 1 and the lower one for Q = 2. Using the 
form factor bound state formula (iv) we obtain (up to a constant) 

It is interesting to note that for Re Q\ — > oo we have the relation of order and 
disorder parameter form factors (up to constants) 

lim (0\a 1 (0)\p 1 ,p 2 ,p 3 ,p 4 )™ 11 = (0|At 2 (0)|p 2 ,p 3 ,P4>iii 
which follows from the asymptotic behavior 
F(9 H ) - el " 



The para-fermi field iPq(x): We look for a solution of (i'3)-(v' 3 ) with 



charge 



spin 



Q 

s 



Q 

Q(3-Q)/3 



statistics <r^ <3 



These fields have charge Q = 1,2 and dual charge Q = Q. The only non- 
vanishing form factors with n particles of type 1 are the ones with n = Q mod 3. 
We propose the correspondence of the field and the p-function: 



„V>i 

rnmk 
rnmk 



pexp E^i-^E^ 



i=i 



i=i 



with < 



pexp [J2 zi + E Ui 
\i=i i=i 



ra = 3/ + Q 
m = Z + 1 

A; = / + g-i 
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where p = y/u Q( * Q Again the normalization constants N n follow from (29). 

Examples for Q = 1,2: The form factors of the para-fermi field ipi{x) for 1 
and 4 particles of type 1 are 

Ft {9) = <0| Vi(0) \ P )i = NJ^iO^) = el e 
i*ii(fi) = (0|^i(0)|pi,p 2 ,P3,P4>Sii = W4 3 i(fi,P* 1 ) II 

l<j<j<4 

4 

= const <f 5 £ * V e 0i +^ TT F ^ 

l<i<j<4 Sillh I (% - ^) Sillh I (% + ^) 

and the one of the para-fermi field ip2(x) for 2 particles of type 1 is 

F$(6) = (0| V 2 (0) |pi,p 2 )n = N 2 I 211 (9,p^)F(0 12 ) 

= const e^ 1+f?2 ) F( ^ 12) . 

sinh |(0i2 - irj) sinh \{6 12 + irj) 

All these examples agree with the results of [18]. 

The higher currents J^(x): We look for a solution of (i'3)-(v' 3 ) with 

charge Q =0 
spin s = L ± 1 
statistics cr^ = 1 . 

Since the currents are Z(3)-charge-less the only no n- vanishing form factors 
with n particles of type 1 are the ones with n = 0mod3. We propose the 
correspondence of the currents and the p-functions for L G Z 



\i=l / \i=l i=l 




Note that for this case the function p in (28) is non- vanishing, however, it does 
not contribute because I nm ,m(9, 1) = for n — 3m. The proof of this fact is 
similar to that one in appendix B. The higher charges Q L — J dxJ^(x) satisfy 
the eigenvalue equations 



n 



QL-T, eLei ) \Pl,.-.,Pn) m = 0. 



i=l 



Obviously, from J±i(x) we obtain the energy momentum tensor. 
Examples: The form factors of energy momentum tensor that is Jf;{x) for 
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L — il for and 3 particles of type 1 are 

F j l = (0| J^(0)|0) =0 

F J A{0) = (0\J£(0) \pi,P2,Ps)tu = csnhn&p^) U 

i<«<i<3 

= ±const (e ±ei + e ±d * + e ±$s ) (e Lei + e Lf?2 + e Le *) 
x TT 

l<*<i<3 Sillh fflii - iT i) Sillh fflij + ^) ' 

By (iv) we obtain the bound state form factor (up to a normalization) for 
L = ±1 

^(0) = icons* e3 (L±1)(< ' 1+ * a) F 1 7 n (0 12 ) . 
Notice that this last expression agrees with the results of [20] when iV = 3. 



4 Z(N)-form factors 



The scaling Z(7V)-Ising model possesses particles of type a = 1, . . . , N — 1 
with Z(iV)-charge Q a = a such that the anti particle of a is a = N — a. The 
bound state fusion rules are (a(3) = a + (3modN, in particular the bound 
state of iV — 1 particles of type 1 is the antiparticle 1. Therefore applying 
iV — 2 times formula (iv) and once (iii) we obtain the recursion relations for 
form factors where only particles of type 1 are involved 

, 9 n ) = 2iF n _ N (9 N+1 , . . . , 9 n ) 

(n \ N-2 

l-af J] S ( Ni) II v^r^ (31) 
i=N+l ) a=l 

where 77 = ^ and the T\^ a = i \Rese = i Va Si a (6)\ 1 ^ 2 are the bound state in- 
tertwiners of the fusion (la) = 1 + a. We will construct the form factors of 
particles of type 1, all the others are then obtained by the bound state formula 
(iv). 

4-1 The general Z(N)-form factor formula 

Following [1] we write the form factors again in the form (1) 

Ffm=KSm n ( 32 ) 

l<i<j<n 



Res . . . Res F n (9 1 , . . . 

N -iN=iri 0i2=iv 
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where minimal form factor function [17] 



dt 2 cosh ±t sinh / , / \\ 

is the solution of Watson's equations (2) with the S- matrix (4). The constant 
is given by (A.l) in appendix A. The K-function K®(0i, . . . , n ) is totally 
symmetric in the rapidities 9 iy 2iri periodic, containing the entire pole structure 
and determines the asymptotic behavior for large values of the rapidities. 
Similar as above we make the Ansatz for the K-functions 

K°(0)=N n I nm (0,p°) (34) 

with the fundamental building blocks of form factors 

(N-l 1 m k 7 \ 

nin / ^f)h(9_,z) P °(0,z) (35) 
k=1 m k l j=1 Jc e _ ti J 

N-l I m k n _ \ 

h(e,z) = n n n - e t ) n - ^ w 

k=l \j=li=l l<i<j<"i k j 

m k mi 

x n n n ^ - %) 

l<k<l<N-l i=l j=l 

where m = (mi, . . . , m^-i) and z = (z k i), k = 1, . . . , N — 1, i = 1, . . . , m k . 
Again the integration contours Cg_ = ^ZCe i enclose the points Oi such that 
R = j Cg dz(p(z — 6). Equations (iii) and (iv), in particular (31) lead to the 
relations 

N-2 N-l 

IJ <f>(z + kiri) J] F(z + kirj) = 1 (37) 

fc=0 k=0 

r(-z)(p (z + in) <p (z) = 1 , x(z)<p(z) = 1 (38) 

/ \ N-l N-2 N-l N-2 N-2 

N n Res0H?) 1] ^( ki v) II F N - h (kir]) = N n _ N 2i II v^S" 

Ve_lr? 1 k=l k=l k=l a=l 

(39) 

where m — 1 — (mi — 1, ... , m^-i — 1). The solution of (37) for is again 

= • u i -\u (40) 

smh ±z smh ±(z + if]) 

if the constant is fixed as in (A.l). The phi-function satisfies again the 'Jost 
function' property (f>(—0)/<f>(0) = S(0). The p-function Pn m (0,z) is analytic 
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in all variables and satisfies: 




Pnm(d,z) is symmetric under 0j <-> 0j 




ff M + 2f 1) fl v .,i)=^(9 1 ,« ; 



'2, 





if # fcfc+1 = ir ] = 2ni/N for A; = 1, . . . , N 



1 



(41) 




= °?Pn-Nrn-Mi^) +P 



where 0' = (6 N +i, ■ ■ ■ ,0 n ), z' = (z ki ), k = 1, . . . , N - 1, % = 2, . . . , m k . In (ii') 
and (iii') af is the statistics factor of the operator O with respect to the 
particle of type 1 and in (v') s is the spin of the operator O. Again p must 
not contribute after integration (in most cases p — 0). 

By means of the off-shell Bethe Ansatz (3) and (35) we have transformed 
the complicated form factor equations (i) - (v) to simple equations for the 
p-functions (i') - (v'). Again one may convince oneself that the form factor 
satisfies (i) and (ii) if h(6_,z_) is symmetric under 0, <-> 9j and periodic with 
respect to 6i — > 0j + 2ni. Not so trivial is again the residue relation (iii) which 
is proved in the following lemma. 

Lemma 2 The form factors given by equations (32) - (36) satisfy the form 
factor equations (i) - (v) if the functions 0, r, x satisfy (37) and (38), the 
normalization constants satisfy (39) and the p-functions satisfy (i') - (v') of 



The proof of this lemma follows the same strategy of the previous Z(3) case. 
Here, however, the essential calculation is much more involved, due to the ex- 
istence more types of particles. Details of this proof can be found in Appendix 



4-2 Examples of fields and p-functions 

We present solutions of the equations for the p-functions (i')-(v') of (41) and 
some explicit examples of the resulting form factors. We identify the fields by 
the properties of their matrix elements. In section 5 we show that the fields 
satisfy the desired commutation rules. This motivates to propose a correspon- 
dence of fields (f>(x) and p-functions p^(9,z). 

The fields ip Q a(x): These fields have the charge Q — 0, . . . , N — 1 and the 



(41). 



c. 
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dual charge Q — 0, . . . , N — 1. We look for a solution of (i')-(v') with 



charge Q mod N 

spin s * = min(Q, Q) - QQ/N (42) 
statistics erf = uo Q 

with a; = e ir? = e 2m ^ N . The phase factor erf is the statistics factor of the 
field iPqq(x) with respect to the particle of type 1. Since the fields carry the 
charge Q the only non-vanishing form factors with n particles of type 1 are 
the ones with n = Q mod N. We propose the correspondence of the field and 
the p-function: 



Q 



^q,q ~ PSS = A>«P I E E zv - % E 0i ] (43) 

\fc=lj=l JV 1=1 



with n = 3/ + Q, I = 0, 1, 2, . . . and 



m fc = / + 1 for k < Q 
m k = I for Q < k 



where p = ^P^~ N+2 ^ n / N \ One easily checks that this p-function satisfies the 
equations (i')-(v') and the requirements (42). The normalization constants N n 
follow from (39). In particular we have for 

Q = the order parameters ctq(x) = i/jqo(x) 
Q = the disorder parameter Pq(x) = i^ q(x) 
Q = Q the para-fermi fields ^q{x) = i[>qq(x) . 



They satisfy space like commutation rules (5), derived in the next section. The 
para-fermi fields iPq{x) are the massive analogs of the para-fermi fields in the 
conformal quantum field theory of [14,15]. One obtains a second set of fields 
iPq q(x) by changing the sign in the exponent of (43). 

The higher currents Jl(x): These fields are charge-less, have bosonic statis- 
tics and spin L ± 1. The only non- vanishing form factors with n particles of 
type 1 are the ones with n = mod N. We propose the correspondence of the 
currents and the p-functions for LgZ 



■ n N—l m 

Jt ^Pn L m = ±J2 e±dl E E eLZkj for n = 3m. 

i=l k=l 3=1 
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The higher charges Ql — J dxJl(x) satisfy again the eigenvalue equations 

(ql-j:^) \Pi,...,p n ) in =o. 



\ 1=1 

T± 



Obviously, from J±i(x) we obtain the energy momentum tensor. 

Examples: Up to normalization constants we calculate for the order param- 
eters (Ti(x) and a 2 (x) 

(O|(Ti(O)|0>i = 1 

F(8 12 ) 



(o\<T2(o)\e 1 ,e 2 yz = 



sinh ±(0 12 - 2ni/N) sinh ±(0 12 + 2ni/N) 
and for the para-fermi fields i/)q(x) and ?/> 2 (:c) 

(0\M0)\0) Q = e £i ^ le (44) 

(i_| r )( 01+e2 ) F ^ , 

(o | v 2 (o) 1^,^)11 = 



sinh §(0 12 - 27ri/iV) sinh ±(0 12 + 2^/^) 

where | 0)q denotes a one-particle state of charge Q and | 6 l 1 ,0 2 )ii a state of 
two particles of charge 1. 



5 Commutation rules 



5. 1 The general formula 



Techniques similar to the ones of this section have been applied for the sim- 
pler case of no bound states and bosonic statistics in [3,27]. A generalization 
for the case of bound states has been discussed in [28]. Here we generalize 
these techniques for the case of more general statistics and also discuss the 
contribution of poles related to the double poles of bound state S-matrices 5 . 
In order to discuss commutation rules of two fields 4>(x) and ip(y) we have 
to use a general crossing formula for form factors which was derived in [4] 
(see also [3]). For quantum field theories with general statistics we introduce 
assumptions on the statistics factor of a field ip(x) and a particle a. It is easy 
to see that for consistency of (ii) and (iii) the condition a^(a)a^(a) = 1 has 
to hold if a is the anti-particle of a. We assume that 

a^(a) = a*{Q a ) (45) 

5 For bound state form factors there are also higher order 'physical poles' (see e.g. 
[29,30,31,32]). 
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depends on the charge of the particle such that ^(Q + Q') = a^(Q)a^(Q'). 
A stronger assumption (which holds for the Z (N)-model) is the existence of 
a 'dual charge' of the fields such that 

a* (a) = u^ Qa (46) 

where |u;| = 1. 

In order to write the following long formulae we introduce short notation: For 
a field 0(x) and for ordered sets of rapidities 9± > ■ ■ ■ > 9 n and 9[ < • ■ ■ < 9' m 
we write the general matrix element of O(0) as 

ot&M := ° ut (Pm(0'J, ■ ■ • I O I a^), . . .,a n (9 n ) ) m (47) 

where 9 a = (9\, . . . , 9 n ) and 9'* = (9[, . . . , 9' m ). The array of indices a = 
(ai,..., a n ) denote a set of particles (a* G {types of particles}) and corre- 
spondingly for (3 (we also write \a\ = n etc.). Similar as for form factors this 
matrix element is given for general order of the rapidities by the symmetry 
property (i) for both the in- and owt-states which takes the general form: 

if # 7 is a permutation of # a and 9^ a permutation of The matrix SaiQ^) 
is defined as the representation of the permutation ^(9^) = # 7 generated 
by the two-particle S-matrices SZH\(9 12 ), for example SZ 3 1 l 2 2 ^ 3 (9 1 ,9 2 ,9 3 ) = 
SZl(9 13 )S£ll(9 23 ) (c.f. [4]). 

We consider an arbitrary matrix element of products of fields O = <p(x)ifj(y) 
and O = ip(y)(f)(x). Inserting a complete set of intermediate states |£y) 7 n we 
obtain 

where = 0(0), ip = -0(0), P a — the total momentum of the state \ 9 JOl ) % a 
etc. and Jg = Ilfcli / Einstein summation convention over all sets 7 is 
assumed. We also define 7! = rL n cJ where n a is the number of particles 
of type a in 7. We apply the general crossing formula (31) of [4] which is 
obtained by taking into account the disconnected terms in (ii) and iterating 
that formula. Strictly speaking, we apply the second version of the crossing 
formula to the matrix element of 

i7r_)C£2l§(g T ,L)S^(£y) (49) 



0f(^i 7 )= E m&Mtsifav*- 
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where p = (p\ p \, ■ ■ ■ , Pij with p = antiparticle of p and &_p — iix^ means that 

all rapidities are taken as & — i (n — e). The matrix 1^(9^ T ,9J) is defined by 
(47) with = 1 the unit operator. The summation is over all decompositions 
of the sets of rapidities &g and 9_ 1 . To the matrix element of ip we apply the 
first version of the crossing formula 

l&jiyM = 4) E s k(LL) WLej c^Via& + sf %) 

(50) 

where we assume that the statistics factor cr^ of the field ip with respect to 
all particles in 7 is the same for all 7 which contribute to (48) (see below). 
Inserting (49) and (50) in (48) we use the product formula 5^(^)5^(0,,, n ) = 
S^{9_ v ,9_ % ). Let us first assume that the sets rapidities in the initial state 9^ 
and the ones of the final state #^ have no common elements which implies 
that also u n 9 a = 0. Then we may use (ii) to get Sff (#„, 9 n ) = 1 and we can 
perform the 6_ u - and ^.-integrations. The remaining ^-integration variables are 
6^ = O^n 9 n , then we may write for the sets of particles q = pu, n = urr and 
7 = purr and similar for rapidities and momenta. Equation (48) simplifies as 



(*Wrt)Mi) 



' 1 10' -a> p 10 1 ■ j S.cj 

x siiQjei^-^-^-^ (51) 



where 



x + i*-,h + in-,h)e~ iPa( - x - y) • (52) 



pT 

The integrand X-^ may be depicted as 



X 



PL 



'(7) 







// 



A 



Similarly, if we apply for the operator product ip(y)(f)(x) again the second 
crossing formula to the matrix element of and the first one the matrix 
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element of ip we obtain equation (51) where Xj^ x is replaced by 

Y S = 4)^^'^ - m ->^ " ^_)C^C^C^ 

xfe(^ + ^_,^,^)e^^ (53) 

which means that only <r^ is replaced by P w by — P w and the integration 
variables Q^hy 9^ — in-. 

No bound states: In this case there are no singularities in the physical strip 
and we may shift in the matrix element of ip(y)(f)(x) (51) with (53) for equal 
times and x 1 < y 1 the integration variables by 0j — > 0j + in- ■ Note that the 
factor e %p ^ x ~y) decreases for < Re 0j < n if x 1 < y 1 . Because P w — > — P w we 
get the matrix element of (j>(x)ip(y) (51) with (52) up to the statistics factors. 
Therefore we conclude 

<j>{x)tj}{y) = ^(y)(f)(x)a^ for x 1 < y 1 (54) 

where = (T ( 7 )/ cr (^)- Using the assumption (45) we have with = I] 7S7 

4) = n = ^(%) = - <w 

7G7 

which is the same for all 7, as assumed above. The last equation follows from 
= — Q 1 and charge conservation which means that the the matrix elements 

in (48) are non-vanishing if Qp + Q<p = Q T Therefore the statistics factor 
of the fields ip with respect to (j) is 

which is in general not symmetric under the exchange of ip and 0. Finally, we 
obtain the space like commutation rules 



l/cr^(Q ) for < y 1 
</>(x)ip(y) = ip(ymx) { (56) 

<jt>(Q f ) for > y 1 



where the second relation is obtained from (54) by exchanging <-> ^ and 
x <r^ y . The same result appears when there are bound states. This is proved 
in appendix D where also the existence of double poles in bound state S- 
matrices is taken into account. 
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5.2 Application to the Z(N)-model 



The statistics factors in this model are of the form (46) a^(a) = u Q ^ Qa where 
is the dual charge of the field ip and Q a is the charge of the particle a, 
therefore = uj~Q*Q<i>. The general equal time commutation rule (56) for 
fields iPqq{x) defined by (43) in section 4 reads as 

! U -RQ — e -2iriRQ/N f Qr x l < 1 
hr 9 win , i i (57) 
U QR = e 2mQR/N for x l > y l 

Notice that in this model we have a more general anyonic statistics. 
Examples: 

(1) The order parameters have bosonic commutation rules with respect to 
each other 

<*Q(x)<TQ>(y) = a Q/ (y)a Q (x) . 

(2) The disorder parameters have again bosonic commutation rules with re- 
spect to each other. 

(3) For the order-disorder parameters we obtain the typical commutation 
rule 



X) 



1 for x 1 < y 1 

U QQ = e 2niQQ/N for x l > y l 

(4) The para-fermi fields have anyonic commutation rules 

M*)My) = My)M*)e e(xl - yl)2mQR/N ■ (58) 

These results prove the commutation rules (5) in the Introduction. 

The 2-point Wightman function: In order to compare these commutation 
rules with the explicit results of the previous section we calculate the 2-point 
Wightman function for the para-fermi fields ipQ and ipN~Q (with spin s = 
Q(N—Q)/N) in 1-particle (charge Q) intermediate state approximation. Using 
the result (44) we obtain 



/riH 
—(o\^)\0)qq(^n- Q (o)\o) + ... 

= ' ^ (My/i(a+-ie)y/i(x--iej)+... 

where v = 2Q(N — Q)/N and x^ = t =F x. This agrees with the commutation 
rule (58), because for t = and x > using the symmetry Q <-> N—Q, x — > — x 
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and translation invariance we obtain 



<0|^(*)^_q(0)|0> 



(01^^)^(0)10) 

e^(0\^ N _ Q (-x)^ Q (0)\0) 

e^(0\^ N _ Q (0)^ Q (x)\0) 



where ((x — it) / (—x — ie)) u = e lnu ^ x ^ 2 has bee used. The asymptotic be- 
havior is obtained from 



for v ^ 0. Therefore the leading short distance behavior is up to constants 
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(0|^q(x)^_q(0)|0>~ (x + -ie) 
(0\i> Q (x)$ N - Q (0)\0)~ (x~-ie) 
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Appendix 

A Some useful formulae 



In this appendix we provide some explicit formulae for the scattering matri- 
ces and the minimal form factors which we frequently employ in the explicit 
computations. The S-matrix of two 'fundamental' particles (i.e. of type 1) is 
[9] 

= = _ exp r * a "X'-w sinh tx 



sinh±(0-^) Jo t sinht 

where is the rapidity difference defined by 

P1P2 = m2 cosh 6 . 

A particle of type a (0 < a < N) is a bound state a = (cci of particles 

of type ctj where a = «i + • • • + ai, in particular a = (1 • • • 1) for all ctj = 1. 

For the scattering of the bound state a and f3 we have [17] 

/-oo dx CQ8hx(l ~ gg) ~ COshx(l - gg) . 

S a/ 3(0) = exp 2 / . 77 , 7T7\ — smhx— . 

Jo x smhxtanh(x/A') «7r 

The minimal form factor functions, which satisfies Watson's equations, are 
obtained from the S-matrix formulae [1] and are given as (f3 > a) 

, dt sinht fl - 4) sinht f-^) / / q\\ 

F^(6)=exv -2 \ N > ,. KN) 1-cosht 1- " 



1 /" c 

a/3 (0) - ex Py o 7* sinh 2 ttanht/iV V' V* i7r , 

in particular[17] 

rfio sinh t (l-i)co S hti f / _ «\\ 



F 1 T(0) = exp/ —2 
jo i 

/■c 

= — i sinh 1 6* exp y 

>smh 5fll r(* + ± + f) r(Jb + i + i-f)Vr(* + i-^ 



sinh t \ \ «7r , 

jtsinhf(l-|) / _ cQsh f /_^ 
i sinh t V \ in 



and with 1 = (N - 1) 

'■oo (it sinh t 2 



in , 



_ ft r (fc + £ + £ + r(fc + | + j - ^) ^ + 1-^) ' 

iir(fc + i-i + ^) r(* + §-±-£) ^r(* + i + i). 
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There are simple relations between the minimal form factors which we essen- 
tially use in our construction which are up to constants 



FT (o + |) {0 - I) oc sinh \ (0 + sinh I ((9 - |) itf»(6>) 

iV-l i, JV-2 / k \ ( k + 1 \ 

II ^iT n (^ + ^2«) oc 1] sinh I f 6> + -2ni J sinh I f + -^2™] 
F^ m (9)F^ m (9 + in) oc sinh ±0 sinh § (0 + 2in/N) . 

In equations (20) and (33) we used the function F(9) = CNF™ n (9) with 
such that the normalizations in (37) and (40) hold. 



B Integrals for the Z(2)-model 

The claim (15) follows from the following lemma 
Lemma 3 For n = 2m + 1 odd and Xi = e° l 

/„te):=/,„,(9,l)-(2 8 )'"- 1 " 2 II ^^T^O- 

Ki<j<n r 



Proof. Again as in the proof of lemma 2 in [33] we apply induction and 
Liouville's theorem. One easily verifies fi(x) = fsix) = 0. As induction as- 
sumptions we take f n ~2 — 0. The functions f n (x) are a meromorphic functions 
in terms of the x^ with at most simple poles at X{ = —Xj since pinchings appear 
for Zk = 9i = 9j ±in. The residues of the poles are proportional to / n _ 2 as fol- 
lows from the recursion relations (iii) for both terms. Furthermore f n (x) — > 
for Xi — > 00. Therefore f n (x) vanishes identically by Liouville's theorem. ■ 

Note that the integrations in the definition (9) of I nm can easily be performed 
and with M — {1, . . . , n} and \K,\—m 

wm)= e n n -^e-- 
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C Proof of the main lemma 



In this appendix we prove the main lemma 2 which provides the general Z(N)- 
form factor formula. 

Proof. Similar as in the proof of lemma 1 we calculate 



mm I N ~ 1 m k r 

Res . • • • Res J^&rfJ = m y m ^ JJ JJ / 
y N -iN=tv e 12 =iv — — mi!...rojv-i! \ r_, ,_ Jc 



mm I N ~ 1 m k r rl-y 

^nm) 

N—l I n m k 



mi ! . . . m N -i ! \ ^ Jc eJ R 



x n n n<t>(zkj-8i) n r( Zki - Zkj ) 

k=l \i=N+lj=2 2<i<j<m k 

m k mi I N—l N m 

x n n n - %) n n n -e i )\r 

l<k<l<N-li=2j=2 \k=li=lj=2 



with r 



/JV-1 \ JV-1 / n _ N 

Res •••Res TT / dz kl ) \[ TT <j>(z k i - 9i) H r(z k i - z kj ) 
e N _ 1N =rr, e 12 =t V \ k=1 J Cl J k=1 \. =1 2< m 



X 



i=2 j=2 



Replacing Cq_ by Cq> where 6' = (9n+i, ■ ■ • , n ) we have used r(0) = r(±ir]) = 
x(0) = x(—irj) = and the fact that the ^-integrations give non-vanishing 
results only for z k \ = 9 k and 9 k+ ±, k = 1, . . . , N — 1. This is because for 
6*i2, • • • , N _ lN — >■ *?7 pinching appears at (z n . . . , , zjv_h) = (0 2 , • • • , #/v) and 
. . . , 0n-i)- Defining the function 



N-l 

f(z u , z N -u) = n n r ( z ki - z kj) n m*** - 

fc=l 2<j<m k l<k<l<N-l 

(mi m k \ 

j=2 i=2 / 



one obtains by means of (38) after some lengthy but straight forward calcula- 
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tion r = 



Res. •••Res / II II " *0 K ■ ■ ■ > ^-11) 



N-2 \ JV_1 /N-l n 



= Res0(-6>) n 0(^) n n 

fc=l / V fe=l i=AT+l 



It has been used that /(..., z .. .) = /(. . . , z, . . . , z — irj . . .) =0 because 

of x(0) = x(—irj) = 0. Using further the defining relation of in terms of F 
(37), the Jost property (27) of the ^-function and the properties (iii') of (41) 
for the p-function we get 



Res . • • • Res I nm (9, V ° nm ) = Res <j>{-0) J[ <p k (kt V ) 

9 N -l N =ir] 0i2=tr] — \0=iri J 



k=l 

N n \ 1 ( n 

Xl n _ Nm _M,Pn-Nrn-l)[U U F(9 H )\ 1 - of ]J S(9 Ni ) 

K k=li=N+l J \ i=N+l 

which together with the relation for the normalization constants (39) proves 
the claim. ■ 



D Proof of the commutation rules 



In this Appendix we prove that we find the same commutation rules for two 
fields 4>{x) and ip{y) when there are bound states poles or even when the 
S-matrix has double poles 6 . 

Bound states: We now show that the same result (54) appears when there 
are bound states 7 which means that there are poles in the physical strip. Let 
7 = (cx/3) be a bound state of a and (3 with fusion angle rp a p which means 
that at 6 a/ 3 = ir? a a the S-matrix S a p(6) has a pole. The momentum and the 
rapidity of the bound state are 

P 1 =P a + Pp 

# 7 = # Q - i(tt - rfi a ) = 9p + i(?r - T]%) 



6 These poles appear typically for bound state - bound state scattering. The case 
of higher order poles may be discussed similarly and will be published elsewhere. 

7 Here we follow the arguments of Quella [28] . 
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where v^ a and 77^ are the fusion angles of the bound states (3 = (701) and 
a = (/97), respectively. 



We start matrix element of ip(y)<j)(x) (given by (51) with (53)). First we con- 
sider the contribution in the sum over the intermediate states where a G u 
and /3 e A. All the particles which are not essential for this discussion will be 
suppressed. Then the function ^ ai g(0 a , nas a P°l e a t 9 a — 9p = iVa/3 sucn 



that by shifting the integration 9 a 
contribution 



9 a + in_ there will be the additional 



- Fies M6 a - i7r_)C a< V a/3 (0 a , 9 p )e 



iP a (x-y) e -iyP /3 



J a — "a 

i 



with 9 a = Op + ^2/3> ^7 = ^/3 + ^(^ ~~ Vpy) an d the fusion intertwiner (see 
e.g. [33]). Next we consider the contribution to the sum over the intermediate 
states where 7GW and (3 = \i. Then the function 0^(0/3, # 7 — in) has a pole 



at 9/3 — 9 1 + in = irf^ such that by shifting the integration # 7 
there will be the additional contribution 



+ in _ 



- Res ■ 

2 7 =6» 7 



# 7 - inSjC^^^y^-^e-™ 1 '? 



= -^MOa ~ in_)V2T%C^(9,y xP «-^ 

with 9 a , # 7 as above and the fusion intertwiner rjg . From the crossing relation 
of the fusion intertwiners 

C<5a-p7 pa /~<77 

1 a/3 - 1 /3 7 ^ 

we conclude that these residue terms form bound state poles cancel. The steps 
of the argument may be depicted as 






J 












Double poles: 8 Form factors have more poles which are not related to 
bound states, they belong to higher poles of the S-matrix. First we consider 
the contribution to the sum over the intermediate states where the particles 
1, 1 e u>, 2 e ji. Again we suppress all particles which are not essential for our 
discussion. Then the function 4> 2 ii(0, — in, 9' — in) has a pole at 9 — 9\ — 

8 This discussion is new. 
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9 + iir/N which correspond to the bound state (21) = 1 with the fusion angle 
7721 = tt(1 - l/N). We shift the integrations 9 — > 9 + m_ and & -> 6*' + m_ 
such that during the shift < Im(# — 9') < e. From the ^-integration there 
will be the additional contribution 

- Res 02ii (9, 9 - in, 9' - i7r)C In Vii(^, 9)e iP ( x -v) e ~ ixP 

2 0=0! 

= - % -(t> u {9 - in/N, 9' - in^r^C 1111 ^^, 9 l )e iP{x ~^ e~ ixP . 

Further the function <pn{9 — in/N, 9' — iir) has a pole at 9' = 9 2 = 9 + — 
3/iV) which correspond to the bound state (11) =2 with the fusion angle 
77^ = n2/N. From the integration there will be the additional contribution 



— ) V2 Res <f> n (9 - iir/N, 9' - i^V^T^C 1111 ^^' , 9 1 )e iP{ - x ~ y) e~ ixP 
2/ e'=e 2 

--(f) 2 (9 - i^/N^T^C 1111 ^^, 9 1 )e i{Pl+P2)(x - y) e- ixP 
2 



This procedure may be depicted as 



1 i x 



Res 




We show this the additional term is cancelled by a contribution to the sum 
over the intermediate states where 2 e u, 2 e A. Then the function i/j 22 (9,9) 
has a pole at 9 — 9 3 — 9 + — 2/N) which correspond to the double pole 
of the S-matrix 



S 22 {9) 



'sin ^ (t + 

2 \VK 



N-2 
N 



2 V VK N 



sin 



TL ( JL 

2 \i-w 



N-2 
N 



sin 



E ( JL 

2 \in 



N-4 
N 



at 9 = Z7r(l — 2/N). From the ^-integration there will be the additional con- 
tribution 



- Res M9 - i7v)C 2 ^- 22 (9, 9)e iP{x - y) e- iyP 



2 0=0 3 

i 



= iM e - iit2/N)C 2 \-i) (i>- n (9 2 , 9 1 )rfC 11 rl 1 ) e iP ^e- iyP . 
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This procedure may be depicted as 




2 - 



Res 
9=e 3 



1 1 11 




The crossing relation of the fusion intertwiners 

2 ^ll 1 2 



r 2 -pllz-ill _ r i22-pll r i_ pill 

1 ill o L — V 1 <) Ulll -i 



11^ 2 

implies that this contribution again cancelled the one above. It has been used 
that the form factor of bound states 22 has a simple pole where the S-matrix 
S 22 has a double pole and the residue is 

Resfe(M) = -i (^i(02,0 3 )lf Cnr* 1 ) . 
8=83 

This may be calculated as follows. By the form factor equation (iv) we have 
Res Res Viin(^i, 2 , 9 3 , 4 ) = 2$j 2 (0 (12) , fl^rfjl* . 

Therefore using the form factor equation (iii) and the definition of the fusion 
intertwiners KesiS = IT we obtain 

8 ^1 2/An ^2(^(12)^(34))rf i r? 1 

= I Res Res Resti Vlll (e 2 ,e 1 ,e±,e 3 )s ri (e 12 )s 11 (e u ) 

l 8 14 =17T #12=1^ 034=^1! 

\ Res Res 2iC n (V>n(0 2 , WitMMM - V>ii(0 2 , 3 )) 



which implies the residue formula used above. This procedure may be depicted 

as 



1> 



6»i 6» 2 6*3 



Res 















1 1\ 


1 



81 82 83 




Res Res 
«12= i '7f 1 e 34= i '?f 1 



81 82 83 




111) |1|1 
2 2T 



Note that the last graph, as an on-shell graph, resembles (half of) the 'box' 
Feynman diagram which was used to investigate the double poles of bound 
state S-matrices (see e.g. [29,30]) 

The general case: Finally we consider the general case that the sets ra- 
pidities in the initial state 9^ and the ones final state 0'p have also common 
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elements. Then after inserting (49) and (50) in (48) there will be S-matrices 
Sff which produce additional poles in the physical strip which would 

produce additional residue contributions while shifting the integration con- 
tours. However, we can remove these S-matrices by using again the crossing 
relation (ii) and move all the lines of common rapidities to the left or right as 
depicted as follows 



(X) 




Then we can apply the procedure as above. 
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